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Abstract 

The Large Hadron Collider has recently discovered a Higgs-like particle having a 
mass around 125 GeV and also indicated that there is an enhancement in the Higgs 
to diphoton decay rate as compared to that in the standard model. We have studied 
implications of these discoveries in the bilinear R-parity violating supersymmetric 
model, whose main motivation is to explain the non-zero masses for neutrinos. 
The R-parity violating parameters in this model are e and b e , and these parameters 
determine the scale of neutrino masses. If the enhancement in the Higgs to diphoton 
decay rate is true, then we have found e •> 0.01 GeV and b t ~ 1 GeV 2 in order to 
be compatible with the neutrino oscillation data. Also, in the above mentioned 
analysis, we can determine the soft masses of sleptons {mi) and CP-odd Higgs 
boson mass (m^). We have estimated that mi ^ 300 GeV and tua ^ 700 GeV. 
We have also commented on the allowed values of e and b e , in case there is no 
enhancement in the Higgs to diphoton decay rate. Finally, we present a model to 
explain the smallness of e and b e . 

PACS numbers: 12.60. Jv, 14.60.Pq, 14.80.Da 
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1 Introduction 



The ATLAS and CMS collaborations of Large Hadron Collider (LHC) have recently dis- 
covered a bosonic particle whose mass being around 125 GeV [lj. The data from the LHC 
is strongly favouring the spin of this bosonic particle to be zero and it is consistent with 
the Higgs boson [2J, which is necessary to achieve the electroweak symmetry breaking. 
The ATLAS and CMS groups have analyzed the decay properties of this Higgs-like parti- 
cle into various standard model fields. An indication for the excess of events in the Higgs 
to diphoton channel as compared to that in the standard model (SM) has been reported. 
Explicitly, by defining the quantity 

o WiPP ^ k ) X BR ( k -> 77)] observ ed m 

[a{pp ^ h)x B R(h -> 77 )] SM ' 1 ) 

where h is the Higgs boson, the ATLAS and CMS had earlier reported that i? 77 = 
1.8 ± 0.5 and 1.56 ± 0.43 [TJ, respectively. The above quoted values for i? 77 have been 
recently updated in March 2013 at the conference Rencontres de Moriond. The ATLAS 
group has claimed i? 77 = 1.65loi3o [3J , which indicates a slight enhancement in the h — > 77 
channel. However, the CMS group has reported that i? 77 could be 0.78^26 or I.II+qJq, 
depending on the type of the analysis [I]. The values quoted by the CMS group imply 
that the discovered Higgs boson is consistent with the SM within the uncertainties. We 
can hope that the future analysis at ATLAS and CMS can resolve the differences in FL fl . 
At this moment, it is worth to analyse by assuming that the discovery made at the LHC 
favours new physics. 

New physics has been motivated by several considerations and some of them are gauge 
hierarchy problem and smallness of neutrino masses. Gauge hierarchy problem can be 
solved by proposing supersymmetry [5] . In supersymmetric models the Higgs boson mass 
can be around the electroweak scale and it is protected from radiative corrections. The 
weakly interacting neutrinos are found to have non-zero masses which should not exceed 
1 eV. The non-zero masses for neutrinos and upper limits on them have been established 
by neutrino oscillation experiments [6] , cosmological observations [7] and /3-decay experi- 
ments [8]. Since the neutrino masses should be smaller than the electroweak scale by at 
least twelve orders of magnitude, the smallness of their masses indicate a new mechanism 
for mass generation. 

To solve both the gauge hierarchy problem and smallness of neutrino masses, bilinear 
R-parity violating supersymmetric (BRPVS) model is a viable option. For a review on 
the BRPVS model, see Ref. [9]. This model is a minimal extension of the minimal super- 
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symmetric standard model (MSSM). In the BRPVS model, additional bilinear terms of 
the forms eLH u and b e LH u are added to the superpotential and scalar potential respec- 
tively. Here L(L) and H U (H U ) are superfields (scalar components) of lepton and up-type 
Higgs doublets respectively. The above mentioned bilinear terms violate lepton number 
and also the R-parity. e is a mass parameter and b e is a mass-square parameter. Provided 
that the e and b e are very small, the masses of neutrinos can be shown to be consistent 
with the observed neutrino oscillation data [101 HH [12] . One may explain the smallness 
of e and b € by proposing additional symmetries [13] or by embedding this model in a high 
scale physics [TJ]. 

The BRPVS model has rich phenomenology [T51 [T6] . In this work we want to study 
the affects of recent discoveries at the LHC on the parameter space of the BRPVS model. 
As mentioned above that the BRPVS model is an extension of MSSM and the additional 
parameters in it are e and b t . Moreover, both e and b e should be very small in order to 
account for the smallness of neutrino masses. As a result of this, in the BRPVS model, 
the contribution to the Higgs boson mass and also to the quantity i? 77 are dominantly 
determined by the MSSM parameters. In order to have light Higgs boson mass nth ~ 125 
GeV, the stop masses should be considerably high as well as large mixing is needed in 
the stop sector [T71 [TBI E] ■ However, parameters in the squark sector do not affect the 
neutrino masses in the BRPVS model. On the other hand, to have _R 77 > 1 it has been 
shown that relatively light stau masses and large left-right mixing in the stau sector are 
required [17] . Essentially, this would mean that the soft parameters of slepton masses 
(mi), higgsino mass parameter (/i) and the ratio of vacuum expectation values (vevs) 
of the two neutral Higgs fields (tan/3) determine i? 77 . We will show that CP-odd Higgs 
boson mass (ttia) also has a role to play in the enhancement of Higgs to diphoton decay 
rate. Shortly below we will explain that the parameters which determine i? 77 can affect 
the neutrino masses in the BRPVS model. It is to remind that in the singlet extension 
of MSSM, enhancement in i? 77 can be made not necessarily with light staus [20] . 

In the BRPVS model, one neutrino state acquires non-zero mass at tree level due to 
mixing between flavor neutrinos and neutralinos [TTJ [12]. The remaining two neutrino 
states acquire masses at 1-loop level due to mixing between sneutrinos and the three 
neutral Higgs bosons pTT|. fT2] . Explicitly, apart from e and b e , the neutrino masses in 
this model are dominantly depended on the neutralino parameters (Mi^, fi, tan/3), mi 
and ttla- From the discussion in the previous paragraph, we can understand that the 
parameters which determine the neutrino masses in the BRPVS model have a role to 
play in the enhancement of Higgs to diphoton decay rate. From this perspective, we can 
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understand that the requirement of i? 77 > 1 can lead to certain allowed values for e and 
b e , which determine the overall scales of neutrino masses. 

Both the ATLAS and CMS groups of the LHC are yet to confirm whether i? 77 > 1 or 
not. Hence we have also analyzed the case _R 77 < 1. In either of these cases we will see 
that the allowed values of e and b t are small, and their smallness can be motivated from a 
high scale physics. While motivating these parameters from a high scale physics, we can 
also predict allowed ranges for m,L, ttla and also about other supersymmetric parameters. 

The paper is organized as follows. In the next section, we give a brief overview of the 
BRPVS model and also describe the neutrino masses in this model. In the same section 
we will also explain the relevant quantities regarding the Higgs boson mass and i? 77 . In 
Sec. 3, we describe our results on e and b e which are compatible with neutrino oscillation 
data and also with i? 77 . We then motivate these results from a high scale physics. We 
conclude in Sec. 4. 



2 The BRPVS model 



The superpotential of the BRPVS model is 



(2) 



where the indices i,j run from 1 to 3. The superfields Q, U and D are doublet, singlet 
up-type and singlet down-type quark fields, respectively. L and E are doublet and singlet 
charged lepton superfields, respectively. H u and H d are up- and down-type Higgs super- 
fields, respectively. As already explained in the previous section, the superpotential terms 
LiH u are the bilinear R-parity violating terms. The addition of these R-parity violating 
terms makes the superpotential of BRPVS model to differ from that of MSSM. Another 
difference between the BRPVS model and the MSSM is that the soft scalar potential in 
the BRPVS model has additional terms which correspond to the LiH u -terms. The form 
of the soft scalar potential in the BRPVS model is 



BRPVS 



soft 



T/ MSSM , 
"soft + 



{b e )iliH u + cc. 



(3) 



Cft SM = l(M 1 BB + M 2 WW + M 3 rg + c.c)+ml u H:H u + m 2 Hd H*H d + 



+{rn 2 Q ) ij Q*Q j + {m 2 u ) ij U*U j + (m 2 D ) l3 D*D 3 + {m 2 L ) ij L*L j + (m 2 E ) l3 E*E 3 
+ UA^ijQiUjHu + (A^ijQiDjHd + (A E ) i3 L i E j H d + b^H u H d + c.c.l . (4) 



The explicit form of the soft terms in the MSSM are given in the form of V r s ^ SM . 
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2.1 Neutrino masses in the BRPVS model 



In this subsection we will describe the neutrino masses, which are generated mainly due 
to the bilinear R-parity violating terms. In fact, these bilinear terms violates lepton 
number, and as a result, the sneutrinos can acquire non-zero vevs. However, without loss 
of generality, we work in a particular basis where the vevs of sneutrinos are kept to be 
zero. 

The e-term of Eq. ([2]) generate mixing between flavor neutrinos (z/,) and higgsino. 
In a basis where ipx 



mixing masses: C 



(B, W 3 , H®, H%, Ui, u 2 , i^3) T , at the tree level we get the following 
^ipJfM N il> N + h.c, where 



M 



N 



M x o m 



m 



(5) 
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Here, g±, g 2 are the gauge couplings corresponding to the gauge groups U(l)y and SU(2)l, 



respectively. The vevs of Higgs scalar fields are defined as: (H®) 



vcosP, {Hi) 



v u = v sin/3, where v = 174 GeV is the electroweak scale. Assuming that 6j are very 
small compared to the TeV scale, at leading order, after integrating away the components 
of neutralinos, we get the neutrino mass matrix as m u = —m T M~ m. However, this 
leading neutrino mass matrix will give only one non-zero mass eigenvalue [TT1 [12] . In 
a realistic scenario we need at least two non-zero mass eigenvalues for neutrinos [Bj. It 
can be shown that the other two neutrino states get non-zero masses due to radiative 
contributions [TR [T2] . At 1-loop level, neutrinos get non-zero masses because of mixing 
between sneutrinos and neutral Higgs boson states [TTJ [12] , and this mixing is driven by 
the soft 6 e -term of Eq. ([3]). It has been shown in Ref. [HI [12] that at 1-loop level, the 
dominant contribution to neutrino masses come from diagrams involving two insertions 
of fe e , provided the tree level mass eigenvalue is dominant over the loop contribution. 
Based on this, below we present the complete expression for neutrino mass matrix. In 
this expression we assume degenerate masses for sneutrinos. 

The neutrino mass matrix in the BRPVS model, up to leading contributions, can be 
written as [TT1[T2] 



{m v )ij = a eiej + ai(b e )i(b e )j, 



(7) 
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where the indices i, j run from 1 to 3. The first term in the above equation is due to 
the tree level effect, which is described in the previous paragraph, and the second term is 
from 1-loop diagrams. The expressions for ao and ai are [HJ [12] 

m 2 7 m^ cos 2 (3 2 . 2 

a = — ^ —-5 ., ., , ■ ™ 7 = cos B w Mi + sin W M 2 , 

fi(m z z rrij sin 2/3 - M 1 M 2 fi) 



a i 



2— (mAToJj {h{m h ,my,m 9 , {m N o)i) cos (a - /3) 

~r COS £5 



1=1 



+h(m H , ma, m 5 , (mjv°)i) sin 2 (a - /?) - J 4 K, m 5 , m p , (m^o)i)) , (8) 

where is the Z boson mass, is the Weinberg angle, and the mjj, m# and itla are 
the light, heavy and pseudo-scalar Higgs boson masses, respectively. The unitary matrix 
Uq diagonalizes the neutralino mass matrix as (Uj M x oU )ij = (m^yo)^-, where {m N o)i 
are the neutralino mass eigenvalues, rrty is the mass of sneutrino field, a is the mixing 
angle in the CP-even Higgs sector. The function I± is given by 

I 4 (m 1 ,m 2 ,m 3 ,m i ) = — =— — s [-^("ii, n^, m 4 ) - 7 3 (m 2 , m 3 , m 4 )], 

mf — m 2 

h(m 1 ,m 2 ,m 3 ) = — - [I 2 (mi, m 3 ) - I 2 (m 2 , m 3 )], 

-i 2 2 

1 mi vn\ 
J 2 (mi,m 2 ) = 77; o 2 2 m — • ( 9 ) 

Since we have assumed degenerate masses for sneutrinos, the neutrino matrix in Eq. 
((7]) generate only two non-zero masses, which is sufficient to explain the solar and atmo- 
spheric neutrino mass scales [6]. By taking the supersymmetric mass parameters to be 
few 100 GeV in ao and a\ of Eq. (J7|), we can estimate the magnitudes of the unknown 
parameters q and (b e )i, in order to have a neutrino mass scale of 0.1 eV. Taking into 
account of the partial cancellations of Higgs boson contributions in a± [T2] of Eq. (J7J), for 
tan/3 ~ 10, we get e« < 10~ 3 GeV and (b e )i ~ 1 GeV 2 . As already described before, the 
estimated magnitudes of e and b € are very small in order to explain the smallness of neu- 
trinos masses, and in this work, we analyze if these estimated magnitudes are compatible 
with the Higgs to diphoton decay rate measured at the LHC. 



2.2 Higgs to diphoton decay in the BRPVS model 

We have already explained before that the BRPVS model is an extension of MSSM, where 
the additional terms are e- and 6 e -terms of Eqs. ([2} and ([3]), respectively. We have argued 
before that both the parameters e and b e need to be very small in order to explain the 
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smallness of neutrino masses. As a result of this, in the BRPVS model, the masses and 
decay widths of Higgs boson states are almost same as that in the MSSM. The leading 
contribution to light Higgs boson mass up to 1-loop level is given by (21] 

ml = ml cos 2 2/3 + — - In — f + - ' * 1 - — ^ , 10 
ft z 4tc 2 v 2 m\ 4n 2 v 2 M 2 \ 12M 2 J v ' 

where m t is the top quark mass, M s = y /mi 1 m i 2 an d X t = A t — \ij tan/3, A t = (Ay) 33. 
Here, 2 are masses of the stops. The second and third terms in the above equation 
arise due to 1-loop corrections from top and stops. The tree level contribution to m^ is ~ 
91 GeV and in order to have m^ ~ 125 GeV, the loop contributions from top and stops 
should be substantially large. As a result of this, the light Higgs boson mass is dominantly 
determined by parameters in the squark sector and the top mass. These parameters do 
not play any role in determining the neutrino masses in the BRPVS model. However, 
to be consistent with the recent Higgs boson mass of ~ 125 GeV, the above mentioned 
parameters should be fixed accordingly in the BRPVS model. It is to remind that the 
loop contribution from top and stop would be maximum if X t = y6Ms- This choice of 
parameter space is known as maximal mixing scenario [21] . In our analysis, which will be 
discussed below, we have considered the maximal mixing scenario in order to have ~ 
125 GeV. 

On the other hand, in order to have enhancement in the Higgs to diphoton decay rate, 
the quantity defined in Eq. (pQ) has a role to play on neutrino masses. Since the dominant 
production for light Higgs boson at the LHC takes place through gluon fusion process, we 
reformulate _R 77 as 

[T(h gg) x BR{h -> 77)] MSSM / ni \ 
77 ~ [T(h gg) x BR(h 7 7)] SM ' 
In the above expression, we have used cr(gg — > h) to be proportional to the decay width 
r(/i — > gg). In the MSSM, the supersymmetric contribution to T(h — > gg) is from squarks, 
while the decay width of h — > 77 gets supersymmetric contribution from squarks, charged 
sleptons, charged Higgs bosons and charginos. For complete expressions, up to leading 
order, for T(h — > gg) and T(h — > 77) in the SM as well as in MSSM, see Ref. [2Tj . 

A scan of parameter space in the MSSM has been done in [T7] and it has been reported 
that to have _R 77 > 1 the masses of staus should be light and the left-right mixing in the 
stau sector should be large. We will show later that i? 77 has some sensitivity to the CP- 
odd Higgs boson mass. The masses and mixing in the stau sector are determined by the 
parameters m 2 L , m 2 R , A%, A* an d tan/3. From the previous subsection, we can notice that 
the magnitudes of e and b € fix the neutrino mass eigenvalues in the BRPVS model. Apart 
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from this, the tree level neutrino masses are depended on the neutralino parameters. Also, 
the 1-loop contribution to neutrino masses are determined by the masses of neutralinos, 
sneutrinos and neutral Higgs bosons. It is to be noticed that the sneutrino masses are 
determined by the soft parameter m 2 L . 

In the previous paragraph we have explained how the neutrino masses in the BRPVS 
model are correlated with the i? 77 . We have studied this correlation and in the next 
section we present our results. Apart from this correlation, one may also study additional 
bounds arising from vacuum stabilization [22], which we leave it for future studies. 

3 Results 

In this section we present our results on the correlation between neutrino masses and i? 77 
in the BRPVS model. We divide this section into three parts. In Sec. 3.1 we describe 
the diagonalization procedure of the neutrino mass matrix of Eq. (J7J, from which we 
obtain expressions for neutrino mass eigenvalues in terms of model parameters. In Sec. 
3.2 we illustrate our method of calculating the i? 77 by varying the model parameters. 
After scanning over model parameters, we can obtain the allowed parameter space of the 
BRPVS model, in order for the neutrino oscillation data to be compatible with i? 77 > 1. 
As stated before that the LHC has not yet confirmed i? 77 > 1, so we make brief comments 
about the possibility of R 11 < 1. From our numerical results we can see that the allowed 
values for e and b e should be very small. We try to motivate the smallness of these values 
from a high scale physics, which we describe in Sec. 3.3. 

3.1 Neutrino mass eigenvalues 

After diagonalizing the neutrino mass matrix of Eq. ([7]), we should obtain mass eigen- 
values as well as the mixing angles. The mixing angles are incorporated in the well 
known Pontecorvo-Maki-Nakagawa-Sakata unitary matrix, t/pMNS> which is the diagonal- 
izing matrix of Eq. ([7]). We parametrize the t/pMNS as it is suggested in [23J. Among the 
three neutrino mixing angles, 612 and 623 are found to be large, whereas, the third mixing 
angle 8 13 is non-zero and is relatively small [23]. From the recent global fit to various 
neutrino oscillation data [23] , we can still choose tri-bimaximal values for # 12 and 623 [2B] . 
Hence, we take sin #12 = sin 623 = ^g. As for the sin #13, at 3cr level, its fitted value 
can be between 0.13 to 0.18 |25j. In our analysis we take sin #13 to be anywhere in this 3a 
range. For simplicity, we assume the Dirac CP-odd phase, 6, and the Majorana phases to 
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be zero. 

From the diagonalization of mass matrix in Eq. (j7J), we obtain the following relation 



m u — ^PMNS m y^PMNS' 



(12) 



where = Diag(m 1; m 2 , m 3 ) and m lj2)3 are the mass eigenvalues of neutrinos. For 
a given set of neutrino mass eigenvalues and mixing angles, the above matrix equation 
can be solved, since it involves 6 relations in terms of 6 unknown parameters (e», (b e )i). 
One possible solution to the above matrix equation is given in [27], in the limit of S13 = 
sin #i3 = 0. Since, now it has been established that S\ 3 7^ [24J, below we describe 
an approximate way of solving the above matrix relation. Although S13 7^ 0, from the 
previous paragraph we can see that S13 ~ 0.1 and hence higher powers of S13 are at least 
one order of magnitude smaller than that of S13. Based on this observation, we can expand 
cos # 13 = a/1 — s?3 ~ 1 — |s?Q H — • . Using this expansion and also fixing the #i 2 and # 23 



to their tri-bimaximal values, we can expand {/pmns in the following way. 
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(13) 



From the above expansion of £/pmns> we can realize that the right hand side of Eq. (JT2]1 
can be expressed as a power series in terms of s\ 3 . Such a matrix relation in Eq. f|12j) can 
be solved if we also express the left hand side of it in a similar power series expansion. 
Hence, we may propose the following series expansions for 6j and (6 e )j. 



Q — Cj,0 + e i,l s 13 + e i,2Stf H ' > 

(6 £ )» = (6 e )i )0 + (&e)i,lSl3 + (&e)i,2Sl 3 + 



(14) 



Here, we can assume that the coefficients e^o, £4,1, Q,2, etc in the expansion for e« have 
same order of magnitude to one another. This also applies to the coefficients in the power 
series expansion for (6 e )j. Below we show one solution for Eq. (I12p . where we solve and 
(b e )i up to 0(s 13 ). 
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Plugging Eq. (jH|) in Eq. f lT2|) . up to (9(si 3 ), we get the following relations. 

[Uqw^Uq].. = 006^06^0 + 01(^)4,0(^,0), (15) 
[{/om^i T + Uxitl^Uq].. = a (e ifi e jA + e^e^o) + ai((6 e )i,o(&e).?,i + (6 e )i,i(6 e )j,o)-(16) 

One solution to the matrix relation in Eq. (TToT) is given below [27] 

ei,o = 0, e 2 , = £3,0 = e, (&e)l,0 = (&e)2,0 = -(&e)3,0 = &e, 

mi = 0, m 2 = 3ai(6 e ) 2 , m 3 = 2a e 2 . (17) 

As already described before, we can understand that m 3 and m 2 are determined by tree 
level and 1-loop level contributions, respectively, to the neutrino masses in the BRPVS 
model. The mass eigenvalue mi has come out be zero, since we have assumed degenerate 
masses for sneutrinos. Using the solution at leading order in Eq. ffTTl) . we can reduce the 
six independent relations of Eq. (fTSj) into five, which are shown below. 

(6 e )i,i = 0, (6 e ) 2 ,i = (& e ) 3 ,i, £2,1 = -£3,1, 

a eei tl + a 1 b e (b e ) 2 ,i, -— -7= = a ee 2 ,i + ai6 £ (6 e ) 2) i (18) 



3V2 U 1,1 3V2 

The last two relations of Eq. (|18p can be solved for infinitesimally many possible values 
of ei,i, £2,1 and (fe e )2,i- We obtain one simple solution by choosing e 2i i = 0. As a result of 
this, for the given set of neutrino mixing angles which we have described above, a solution 
to the matrix relation of Eq. (lT2"j) . solved up to C(si 3 ), is 

6i = e[V25i3 + 0(5 2 13 )], e 2 = e[l + 0(sl 3 )], e 3 = e[l + 0(sl 3 )], (b e ) x = b e [l + 0(s 2 l3 )], 
(b e ) 2 = b e [l - -^s 13 + 0(s 2 13 )], (6 e ) s = b e [-l - -L ai3 + (9(4)]. (19) 

Here, e and 6 e determine the non-zero neutrino mass eigenvalues of the BRPVS model, 
which are given in Eq. ( TT7I) . We believe the procedure described above can be extended 
to solve 6j and (fo e )j up to second and higher order powers of s i3 . 



3.2 Computation of Higgs to diphoton decay rate 

As already explained that the Higgs to diphoton decay rate and the masses of scalar Higgs 
bosons in the BRPVS model are almost same as that in the MSSM. The enhancement 
related to this decay rate, as quantified in Eq. ffTTl) . and also the masses of Higgs bosons 
have been computed with the HDECAY code [28]. In our numerical analysis, we have 
fixed off-diagonal elements of soft mass-squared and A-terms of Eq. (J4j) to be zero, 
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which is also incorporated in the HDECAY code. In order to have the light Higgs boson 
mass to be around 125 GeV, we have fixed (^0)33 = { m u)33 — ( m l>)33 = (800 GeV) 2 , 



(y4[/) 3 3 = y6(mg)33 + /icot/3 and (Ad) 33 = 0. The specific choice for (A(/) 33 has been 
motivated by the maximal mixing in the stop sector [21] . We have fixed the top quark 
mass to be 173.2 GeV. We have not changed the above parameters, since they do not affect 
the neutrino masses in the BRPVS model. Indeed, for the above set of parameters in the 
squark sector, we have almost got rrih ~ 125 GeV, by varying the other parameters in the 
model. As explained before, the neutrino masses in the BRPVS model are determined by 
neutralino parameters (Mi, M 2 , /x, tan/3) and by the masses of neutral Higgs bosons and 
sneutrinos. We have chosen Mi = | tan 2 9\y M 2 and we have varied M 2 from 100 GeV 
to 1 TeV in steps of 100 GeV. As mentioned before that in order to have i? 77 > 1, the 
mixing in the stau sector should be large, which is determined by /i, tan/3 and (^4^)33. In 
our analysis, we have varied fi from 100 GeV to 2 TeV in steps of 50 GeV and tan/3 has 
been varied from 5 to 60 in steps of 5. We have fixed (^4e)33 = 0. As explained before, 
while solving for the neutrino mass eigenvalues, we have assumed degenerate masses for 
sneutrinos. This would imply that (m\)ij = m 2 L 5ij. For right-handed slepton masses, we 
have assumed (m E )ij = m 2 E 5ij and we have fixed m 2 L = m 2 E . We vary the parameters 
and ttla, which determine the sneutrino mass as well as the masses of CP-odd and heavy 
Higgs bosons. 

In the previous paragraph we have specified parameters of the model in order to 
compute rrih and FL fl . In fact, for fixed values of and m^, we scan over the neutralino 
parameters. In the scanning procedure, we have demanded that some constraints need 
to be satisfied. Among these, we have applied constraint from the muon anomalous 
magnetic moment, (g — 2) M [29|. The current world average value of (g — 2) M differs from 
its corresponding SM value by about 3a [29]. This discrepancy in (g — 2) M is quantified by 
Aa^. In the MSSM, at 1-loop level, Aa M gets contribution from neutralino-charged slepton 
and chargino-sneutrino loops [30] . Since we have justified before that in the BRPVS model 
the additional parameters q and (b e )i are very small, hence the contribution to Aa M in 
the BRPVS model is almost same as that in the MSSM. As a result of this, we have used 
the above mentioned loop contributions to Aa M [30] in our numerical analysis. Below we 
describe the four constraints which we have applied in our scanning procedure. 

(i) m h should be in the range of 123 to 127 GeV, 

(ii) either i? 77 > 1 or i? 77 < 1, 

(iii) masses of sleptons should be greater than 100 GeV, 
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(iv) Aa M should be in the range of (1.1 — 4.7) x 10 9 . 

The constraints (i), (iii) and (iv) have been applied in every case. Regarding the constraint 
(ii), we will specifically mention below whether i? 77 > 1 or i? 77 < 1 has been applied. For 
those points in the parameter space which satisfy the above four constraints, we calculate 
e and b e which determine the neutrino masses through Eq. (fTTj) . We have chosen the 
following values for neutrino mass eigenvalues in order to be consistent with the neutrino 
oscillation data [6]: 



Here, the solar and atmospheric mass scales (central values), from a global fit to neutrino 
oscillation data [25]) are given as Am s 2 ol = 7.62 x 1(T 5 eV 2 and Am^ = 2.55 x 1(T 3 eV 2 , 
respectively. 

In Fig. 1 we have shown allowed values of e and b e for different values of and rriA- 
In Fig. 1(d) we have applied the constraint i? 77 < 1, while in other plots of Fig. 1 the 
constraint i? 77 > 1 has been applied. In Fig. 1(a) has been kept to a very low value of 
200 GeV, and in Figs. 1(b) and 1(c) = 300 and 400 GeV, respectively. As explained 
earlier, the allowed points in these plots are satisfied by the requirement _R 77 > 1. We 
have found that for = 150 GeV, the constraint _R 77 > 1 is not satisfied. From the 
plots of Fig. 1(a) — (c), we can notice that the most likely value of e is > 0.01 GeV. This 
value of e is at least one order larger than its expected value from the neutrino masses, 
which is described at the end of Sec. 2.1. For i? 77 > 1, the lowest value of e can be found 
in Fig. 1(a), which is ~ 0.007 GeV, and at these points M 2 should be as low as 100 GeV. 
On the other hand, in future, if LHC has not found any excess in the Higgs to diphoton 
decay rate, then from Fig. 4(d) we can notice that we can satisfy the neutrino oscillation 
data for e between about 10~ 4 to 0.1 GeV. In Fig. 4(d) we have fixed mi = 400 GeV. By 
decreasing mi, the allowed space for e and b e is slightly different from that of Fig. 4(d). 
Hence, from the measurement of Higgs to diphoton decay rate, if i? 77 > 1 then e should 
be at least ~ 10~ 2 GeV. Otherwise, if i? 77 < 1 then e can be as low as ~ 10~ 4 GeV. 

From the plots of Fig. 1(a) — (c), we can notice that when we increase by keeping 
rriA fixed, the average value of e is increasing. By increasing m^, the lower limit on tan/3 
and /i will increase in the case of i? 77 > 1, which we will describe below. As a result of 
this, the quantity ao of Eq. (JHJ), which is inversely related to e, will decrease. 

In Sec. 2.1, from the neutrino mass scale we have estimated that b e is ~ 1 GeV 2 . In 
the case of i? 77 > 1, from the plots of Fig. 1(a) — (c), we can notice that for > 700 
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Figure 1: Allowed values of e and b e . In plots (a), (b) and (c) we have fixed tul — 200, 
300 and 400 GeV, respectively, and we have applied _R 77 > 1. In plot (d), i? 77 < 1 has 
been applied and mi = 400 GeV. In each of these plots red, green and blue color points 
represent ttia = 1000, 700 and 500 GeV, respectively. 



GeV, b e can be in the range of 0.5 to 2 GeV 2 . For i? 77 > 1, the lowest value of b t has 
been found to be about 0.05 GeV 2 , which can be seen in Fig. 1(c). In Fig. 1(d) we have 
i? 77 < 1, and the allowed value of b e can be around 1 GeV 2 by appropriately choosing the 
uia- For instance, if we have to achieve e ~ 10~ 3 GeV and b t ~ 1 GeV 2 , then uia should 
be < 500 GeV. 

In Figs. 1(a) — (c), for a given value of tul, b e is increasing with tua- The reason for 
this is that b e is inversely related to a±, and from Eq. (jHJ) we can understand that a\ 
decreases with uia- Similarly, from Figs. 1(a) — (c), by keeping itla fixed, we can notice 
that the average value of b e is decreasing with m^. We will shortly explain below that by 
increasing mt the lower limit on tan/3 and fi will increase in the case of i? 77 > 1. Hence, 
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although the function J 4 of a\ decreases with increasing m^, the factor eos 1 2 -5 in a\ will 
compensate this decrease, and the net result is that a\ increases with m^. 

In Fig. 2 we have plotted allowed values of fi and tan 0. The points in Fig. 2(a) are 
allowed by the constraint i? 77 > 1, whereas the points in Fig. 2(b) satisfy _R 77 < 1. In 
the plots of Fig. 2 there are no allowed points for tan (3 = 5. We have found that for 
such a low tan/3 the mass of light Higgs boson is below 123 GeV and hence do not satisfy 
the constraint (i). In Fig. 2(a), we can notice that by increasing the value of mi, the 
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Figure 2: Allowed values of /1 and tan/3 for tua = 1 TeV. In plots (a) and (b) we have 
applied the constraints i? 77 > 1 and i? 77 < 1, respectively. The red, green and blue points 
are for m L = 400, 300 and 200 GeV, respectively. 

lower limit on fi and tan /3 would increase. We may understand this from the fact that the 
stau masses should be as light as possible, and moreover, the mixing in the stau sector 
should be large in order to have i? 77 > 1 [T7j. Hence, by increasing the soft mass m L , the 
quantity /i x tan/3 should proportionately be increased in order to decrease the lightest 
stau mass and also to increase the mixing in the stau sector. From Fig. 2(a), we can 
notice that for a specific value of tan /? the allowed value of fi lies in a certain range. We 
have seen that the lower and upper limits of this range are restricted by the constraints 
(ii) and (iii). For instance, for = 300 GeV and tan/3 = 30, the allowed range for fi is 
from 1050 to 1500 GeV. In this case, for fi < 1050 GeV we may not satisfy i? 77 > 1. On 
the other hand, for \i > 1500 GeV the lightest stau mass becomes less than 100 GeV. In 
Fig. 2(a) we have fixed = 1 TeV. By decreasing m^, we have found that tan/3 is not 
restricted, however, for each tan/3 the corresponding lower limit on /i will increase. To 
illustrate this point, by considering the case of ttia = 700 GeV, = 300 GeV and tan/3 
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= 30, the allowed range for /i has been found to be between 1150 to 1500 GeV. Hence, 
these results indicate that by decreasing rriA the i? 77 value will decrease. 

As stated before, in Fig. 2(b) we have applied the constraint i? 77 < 1. In this plot we 
can see that fi can be as low as 100 GeV. Numerically, we have noticed that i? 77 increases 
with /i and hence after a certain large value of fi, i? 77 < 1 may not be satisfied. In the 
case of mi = 400 GeV, in Fig. 2(b), for tan/3 = 10 and 15, large values of /i are not 
allowed by the constraint (iv). In fact, allowed points in Fig. 2(b) indicate that i? 77 < 1 
can be satisfied for large tan /3 and relatively large \i values. For these large values of fi 
and tan/3, the calculated values of e can be as high as 0.1 GeV, which can be seen in Fig. 
1(d). For low tan/3 and moderate values of /i, e can be < 10~ 3 GeV. 

In Fig. 3, we show the correlation between enhancement in the Higgs to diphoton 
decay rate (-R 77 ) and the bilinear parameter e. From Figs. 3(a) — (c), we can observe that 
for a low value of = 200 GeV, the maximum value for i? 77 is ~ 1.1. As noted before, 
in the case of _R 77 > 1, the lowest value of e is ~ 0.007 GeV, which is found for m L = 
200 GeV. For this lowest value of e the i? 77 value is barely greater than 1.0. From Figs. 
3(a) — (c), we can notice that the maximum value of i? 77 is increasing with mi- We have 
stated before that by increasing m^, the values of fi and tan/3 would increase. For large 
values of \i and tan/3, the coupling strengths of staus to the light Higgs boson would 
increase [17J. As a result of this, _R 77 is increasing with mi- The maximum values of i? 77 
in Figs. 3(b) and 3(c) are 1.41 and 1.92 respectively. We may increase _R 77 to more than 
2.0 by increasing to 500 GeV. But in order to have large mixing and lower masses in 
the stau sector, we have to proportionately increase [i and tan /3. In this work we have 
scanned \i and tan /3 up to 2 TeV and 60, respectively, and have not considered cases of 
■mi > 500 GeV. However, from the above mentioned arguments, for i? 77 > 1, we can 
speculate that by increasing to 500 GeV the value of e would be around 0.1 GeV. 

In Fig. 3(d), we have applied the constraint i? 77 < 1. We can notice from this plot 
that for e ~ 10~ 3 GeV, _R 77 is different for different values of tua- From this perspective 
we can argue that, if i? 77 < 1 is found to be true, then a precise measurement of i? 77 can 
be used to determine e and m^. 

3.3 Smallness of e and b e 

In this subsection, we try to motivate the smallness of e and b e from a high scale physics. 
Essentially we will explore what the Higgs to diphoton decay rate can tell us about the high 
scale physics parameters. As it is noted in [T4], by assuming supersymmetry breaking at 
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Figure 3: Allowed values of e versus i? 77 . In the plots (a), (b) and (c), the value of m L 
= 200, 300 and 400 GeV, respectively. In these plots we have applied the constraint 
_R 77 > 1. In plot (d), mi = 400 GeV and the constraint i? 77 < 1 has been applied. In 
each of these plots red, green and blue color points represent ttia = 1000, 700 and 500 
GeV, respectively. 
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an intermediate energy scale A ~ 10 11 GeV, we can explain the //-parameter, soft terms in 
scalar potential as well as e- and b e - terms. Here we briefly describe important ingredients 
from Ref. [Hj. By introducing SM gauge singlet superfields S, X± and X 2 , we may write 
the superpotential as 



w = a2s + Th***' + jL"""- + M LH » + " • • (21) 



Here, Mp = 2.4 x 10 18 GeV is the Planck scale. There can be 0(1) couplings in the 
above terms, which we have neglected. In the above equation we have written only the 
necessary terms for our purpose, and these terms can be justified by introducing additional 
symmetries, say gauged U(l)'. S must be singlet under this additional U(l)', but X 12 
can be charged under U(l)'. The vevs of the scalar components of these SM gauge singlet 
superfields can be arranged as [14] : (S) ~ Mp, (Xig) = Ai )2 ~ A. The first term of Eq. 
f[2Tj) breaks supersymmetry spontaneously by acquiring an auxiliary vev to S which is of 
the order of A 2 . This axillary vev can generate soft terms in the scalar potential with mass 
parameters m so f t ~ Here, the generation of soft terms in the scalar potential is based 
on the Polonyi mechanism [31]. The scalar vevs of X\ i2 generate the \x- and e-parameters 

A 2 A 3 

which are // ~ and e ~ j^. Here, we have followed the Kim-Nilles mechanism for 
the generation of //-term [32] • The auxiliary vevs of X 12 can generate and b e which are 

, AfA 2 , , A?A| 

fr M ~ isr- and b t ~ -jg. 

In the previous paragraph, we have given motivation for the generation of e and b e 
parameters as well as other supersymmetric parameters from a high scale physics. Now 
we have to fix the high scale physics parameters in order to fit the low energy data. Since 
we expect m soft ~ /i ~ TeV, for A ~ A 1 ~ 0.5 x 10 11 GeV we can explain the TeV scale 
masses for supersymmetric fields. If we take A 2 ~ 10 11 GeV we can get e ~ 1CT 3 GeV. 
From Figs. 1 and 3, we can say that a value of e ~ 10~ 3 GeV is consistent with i? 77 < 1. 
In order to achieve i? 77 > 1, e should be > 0.01 GeV. Hence, by taking A 2 ~ 3.9 x 10 11 
GeV we can get e ~ 10~ 2 GeV. So, if i? 77 > 1 then there is a little hierarchy between A x 
and A 2 , otherwise, this hierarchy can be reduced. 

In future, if LHC has found that i? 77 is significantly larger than 1.0, then from Figs. 
3(a) — (c) we can say that should be larger than about 300 GeV. For tjil between 
300 to 400 GeV, from Figs. 1(b) and 1(c) we can have b e ~ 1 GeV 2 , provided uia is ^ 
700 GeV. Now, for Ai ~ 0.5 x 10 11 GeV and A 2 ~ 3.9 x 10 11 GeV, we can get 6 M ~ 
8 TeV 2 and b e ~ 1 GeV 2 . Hence, for the case of i? 77 > 1, we can motivate consistent 
supersymmetric mass spectrum and 0.1 eV scale for neutrino masses from a high scale 



17 



physics by proposing two different intermediate scales. The hierarchy between these two 
scales should be at least 8. 

If there is no enhancement in the Higgs to diphoton decay rate, from Fig. 1(d), we can 
notice that e can be between ~ 10~ 3 to 0.1 GeV. From the above discussion, to achieve e ~ 
0.1 GeV from high scale physics, there need to be little hierarchy between the intermediate 
scales A 12 . This hierarchy can be minimal for e ~ 10~ 3 GeV. For e ~ 10~ 3 GeV, in Fig. 
1(d), b e can be around 1 GeV 2 if rriA ~ 500 GeV. For this set of values, from Fig. 3(d), 
we can notice that _R 77 is little less than 0.9. Hence, if we believe in the motivation of e 
and b t from high scale physics, the high energy scales in this scenario depend on the value 
of _R 77 . Moreover, in the above described analysis, we can also estimate and rriA by 
knowing the i? 77 . So the future runs at the LHC can give us clue about this scenario by 
precisely finding the R yi . 

We make short comments about measuring the parameters e and b t in experiments. 
Since both these parameters indicate that R-parity is violated, a consequence of that is 
that the lightest supersymmetric particle (LSP) is unstable. Depending on the parameter 
space, either the lightest neutralino or the lightest charged slepton can be the LSP in this 
model [16] . The decay life time of LSP is determined by e and b e . Hence the signals of 
the decay of LSP in this model should give an indication about these bilinear parameters 
[16], from which we can verify the neutrino mass mechanism and also its correlation to 
the Higgs to diphoton decay rate. 

4 Conclusions 

Recently the LHC has discovered a Higgs-like particle whose mass is around 125 GeV. 
It has also been indicated that there is an enhancement in the Higgs to diphoton decay 
rate as compared to that in the SM. We have studied implications of these discoveries in 
the BRPVS model. This model is a minimal extension of the MSSM where the bilinear 
terms eLH u and b t LH u are added to the superpotential and scalar potential of the model, 
respectively. The main objective of this model is to explain the smallness of neutrino 
masses, where the neutrino mass eigenvalues can be shown to be dependent on neutralino 
parameters, soft mass for charged sleptons (mi) and CP-odd Higgs boson mass (ttia) 
[TT| [12], apart from the bilinear parameters e and b t . 

In our analysis, we have scanned over the neutralino parameters and varied and rriA 
accordingly. We have also fixed the soft masses for third generation squarks, in order to 
have light Higgs boson mass to be around 125 GeV. We then have studied implications of 
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enhancement in the Higgs to diphoton decay rate (-R 77 ) in the BRPVS model. Explicitly 
we have found that in order to be compatible with i? 77 > 1 and the neutrino oscillation 
data, the unknown bilinear parameter should be e > 10~ 2 GeV. We have also found that 
to achieve _R 77 between about 1.5 to 2.0, should be between 300 to 400 GeV, provided 
/i and tan /3 are scanned up to 2 TeV and 60 respectively. We have not obtained specific 
bounds on m^. However, from the order of estimations we expect b e to be around 1 GeV 2 
and to achieve this with the above mentioned mi, rriA can be in the range of 700 to 1000 
GeV. 

Since i? 77 > 1 is not yet confirmed by LHC, we have also analyzed the case of _R 77 < 
1. In this later case, we have found that e can be between ~ 10~ 3 to 0.1 GeV. The 
corresponding b e can be around 1 GeV 2 by appropriately choosing wia to be from 500 to 
1000 GeV. Moreover, we have also found that _R 77 can be as low as 0.85. 

From the above discussion, we can notice that e and b e need to be very small in GeV 
units. We have motivated smallness of these two parameters from a high scale physics, 
and at the same time we have also explained the TeV scale masses for supersymmetric 
fields. We have found that to explain e ~ 10~ 2 GeV and b e ~ 1 GeV 2 there need to be 
two different intermediate scales (~ 10 11 GeV) with a hierarchy of factor of 8 between 
them. 
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